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Abstract 

The dynamics of solitons of the nonlinear Schrodinger equation under the 
influence of dissipative and dispersive perturbations is investigated. In partic¬ 
ular a coupling to a long-wave mode is considered using extended Ginzburg- 
Landau equations. The study is motivated by the experimental observation 
of localized wave trains (‘pulses’) in binary-liquid convection. These pulses 
have been found to drift exceedingly slowly. The perturbation analysis re¬ 
veals two distinct mechanisms which can lead to a ‘trapping’ of the pulses 
by the long-wave concentration mode. They are given by the effect of the 
concentration mode on the local growth rate and on the frequency of the 
wave. The latter, dispersive mechanism has not been recognized previously, 
despite the fact that it dominates over the dissipative contribution within the 
perturbation theory. A second unexpected result is that the pulse can be 
accelerated by the concentration mode despite the reduced growth rate ahead 
of the pulse. The dependence of the pulse velocity on the Rayleigh number 
is discussed, and the hysteretic ‘trapping’ transitions suggested by the per¬ 
turbation theory are confirmed by numerical simulations, which also reveal 
oscillatory behavior of the pulse velocity in the vicinity of the transition. The 
derivation and reconstitution of the extended Ginzburg-Landau equations is 
discussed in detail. 
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I. INTRODUCTION 


The appearance of stable localized structures in spatially homogeneous extended sys¬ 
tems constitutes an exciting feature of various pattern-forming systems. Although the lin¬ 
ear instabilities arise from extended modes, various nonlinear mechanisms can lead to the 
confinement of the structure to a small part of the system, the size of which is independent 
of the system size. 

In systems exhibiting steady structures the competition between structures with two 
or more different wave numbers can lead to stable domains of one wave number within a 
background of a different wave number 0-^. Such structures have been observed in con¬ 
vection experiments in narrow channels |ll]. They are expected to be particularly accessible 
in parametrically driven surface waves in ferrofluids |^. The stability of domain structures 
is due to an oscillatory interaction between the domain walls separating the domains 
Zig-zag structures constitute a two-dimensional analog of these domains p . 

A different class of localized structures has been found in systems exhibiting a parity¬ 
breaking bifurcation from a steady structure to traveling waves MM- Instead of extended 
traveling waves one observes that the parity (reflection symmetry) of the structure is broken 
only over a stretch of a few wavelengths. Within that part, which is drifting through the 
systems, the structure is traveling, whereas the remaining structure is motionless. Here the 
localization is connected with the fact that generically the extended traveling waves are 
unstable to side-band perturbations right at their onset [|Tl| -[l^. Remarkably, the parity¬ 
breaking bifurcation need not be subcritical for this localization to occur |Tl| , [l^ . 

In this paper localized traveling waves are investigated which arise in the convection in 
binary liquid mixtures. This system has been the subject of extensive experimental and 
theoretical studies. In quasi-one-dimensional set-ups ‘pulses’ of traveling waves have been 
observed with lengths ranging from a few up to 15 wavelengths ||T^-[l8|. In sufficiently 


homogeneous systems they drift extremely slowly through the system [|^. Outside the 
pulses the fluid is in the conductive, motionless state. In two-dimensional systems such 
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pulses have also been observed but were always found to be (very weakly) unstable [jT^. A 
more detailed discussion can be found in a recent review m and in pl| , |22[| . 

Theoretically, the pulses in binary-mixture convection have posed various intriguing prob¬ 
lems. These are related to the mechanism of localization, the extremely slow drift of the 
pulses and their stability. Numerical simulations of the Navier-Stokes equations have pro¬ 
vided quite some insight into the physics of the localized (and the extended) waves and in 
particular into the role of the concentration held pl|-^. In order to get a detailed un¬ 
derstanding of the mechanisms leading to the localization of the traveling waves and the 
slowing-down of the resulting pulse it is highly desirable to study the problem also analyt¬ 
ically. However, in particular due to the subcritical nature of the bifurcation to traveling 
waves, no rigorous asymptotic reduction of the Navier-Stokes equation to a set of simpler 
evolution equations describing pulses appears available. Nevertheless, asymptotic analysis 
allows the derivation of equations of the Ginzburg-Landau type which capture and elucidate 
the essential mechanisms of the system. 

Various authors have shown that the dispersive character of the waves can lead to the 
localization of the waves. In the limit of strong dispersion the complex Ginzburg-Landau 
equation describing the backward Hopf bifurcation to traveling waves was considered as 
a dissipatively perturbed nonlinear Schrodinger equation It was shown that the 

dissipative perturbations could select a stable solitary wave from the continuum of soliton 
solutions. In the opposite limit of weak dispersion the interaction of fronts connecting the 
conductive and the convective state were considered. It was shown that dispersion can lead 
to a spatially oscillatory contribution to the interaction, thus allowing stably bound pairs of 
fronts forming a solitary wave p8|-^. 

In both analytical approaches serious qualitative disagreements with the experimental 
observations remain. First, the theoretical drift velocity of the pulses p2|-p^ is found to 
be by a factor of 20 to 40 and more larger than that observed in experiment ||T8| , |35| . In 
particular, the theoretical pulse velocity is always in the same direction as the phase velocity 
of the waves. Experimentally, however, backward traveling pulses are found as well [^. It 
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has been argued that this discrepancy could be eliminated if nonlinear gradient terms were 
kept in the Ginzburg-Landau equation 1^^. This would, however, require that the linear 
group velocity be compensated for by higher-order correction terms, which would constitute 
a quite delicate balance. Moreover, within the Ginzburg-Landau equation small and large 
pulses coexist with the large pulse being stable and the small one always being unstable 
P^,|30[|. In recent experiments, however, only a single (stable) pulse was found for strongly 
negative separation ratio 'ip and two pulses for weakly negative V’- In the latter case it was 
the shorter pulse which was stable whereas the longer pulse was unstable . 

As pointed out in previous papers mass diffusion in liquid mixtures is extremely slow as 


compared to heat diffusion This introduces an additional slow time scale connected to 
a concentration mode and leads to the break-down of the complex Ginzburg-Landau equa¬ 
tions already for quite small amplitudes of convection. Based on this observation extended 
Ginzburg-Landau equations were introduced which capture this aspect p7|,^. Numerical 
simulations of these equations showed that the additional mode can lead to a consider¬ 
able slow-down of the pulses [^|. An analysis of the interaction of fronts showed that 
this mode introduces a new mechanism for localization, which is independent of dispersion. 
Strikingly, it depends strongly on the direction of propagation of the pulse [Q. Numeri¬ 
cal simulations of the extended equations revealed that more than two pulse solutions can 
coexist. In that case the longer pulse can be unstable in agreement with the experimental 


results 41 


In the present paper I investigate the extended Ginzburg-Landau equations analytically 
in detail. I consider the strongly dispersive case and study in particular the effect of the 
concentration mode on the perturbed soliton solution. In certain limits the concentration 
mode can be eliminated and leads to a change of the coefficients in the Ginzburg-Landau 
equation. In particular, it introduces nonlinear gradient terms. Most notably, its contri¬ 
bution to the coefficients depends on the velocity of the pulse. This leads to a strongly 
nonlinear dependence of the pulse velocity on the parameters. The calculation shows that 
there are two different mechanisms that can lead to a slow-down (‘trapping’) of the pulse. 
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As found in numerical simulations of the Navier-Stokes equations pT[ the concentration 
mode strongly reduces the buoyancy of the liquid ahead of the pulse. This can - but need 
not - slow down the pulse. There is in fact the possibility that the pulse is sped-up by the 
concentration mode. This is due to the fact that the variation in buoyancy also leads to a 
change in the wave number of the pulse which is directly related to its velocity. The second 
mechanism does not rely on a change in the buoyancy. It arises from the change in frequency 
of the wave due to the concentration mode and can also induce ‘trapping’. This was not 
anticipated from the numerical simulations . 

This paper is organized as follows. In sec.|| the extended Ginzburg-Landau equations 
are presented. Their validity and their relevance to other systems exhibiting an interaction 
between short-wave and long-wave modes is discussed in some detail. The concentration 
held generated by a pulse within this framework is described in sec. B The main results 
for the inhuence of the concentration mode on the pulse and its velocity are obtained in 
sec. using perturbation theory around the soliton of the nonlinear Schrodinger equa¬ 
tion. Numerical simulations supporting the conclusions based on the perturbation result are 
presented in sec.|V[ In two appendices a slightly more general derivation of the extended 
Ginzburg-Landau equation from the Navier-Stokes equations is presented, which leads to 
quantitative agreement with the usual Ginzburg-Landau equation in a suitable limit. Their 
self-consistency is conhrmed using Ward identities. 


II. THE EXTENDED GINZBURG-LANDAU EQUATIONS 


As indicated in the introduction, the pulses of the complex Ginzburg-Landau equation 
differ in various aspects qualitatively from those observed experimentally. The numerical 
simulations of the Navier-Stokes equations indicate that certain modes of the concentration 
held are independent of the convective amplitude |^. In the complex Ginzburg-Landau 
equation the concentration held is, however, adiabatically slaved to the convective amplitude. 
The concentration held can become an independent dynamic quantity only if its time scales 
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are slow, i.e. of the same order as those of the convective amplitude. In [^] it was recognized 
that in binary-liquid convection this is indeed the case since the mass diffusion is by a factor 
100 slower than heat diffusion. Since the Ginzburg-Landau equation is only valid as long as 
the time scales of the convective amplitude are slow as compared to all other time scales, 
the slow mass diffusion implies that the regime of validity of the Ginzburg-Landau equation 
is limited to very small amplitudes. 

To capture the additional dynamics of the concentration held I have considered the 
Navier-Stokes equations in an expansion in which the Lewis number, which characterizes 
the ratio of mass diffusion to heat diffusion, is of the same order as the growth rate of the 
convective amplitude 1^^. In the relevant limit of vanishing mass diffusion two types of 
critical modes arise: the usual convective mode A, which involves the stream function, the 
temperature and the concentration held, and an inhnite number of critical modes which are 
all associated with the concentration held alone: any concentration mode which is indepen¬ 
dent of the horizontal coordinate has vanishing growth rate. Strictly speaking, this implies 
that in a weakly nonlinear theory this inhnite number of modes would have to be kept. 
However, since all these modes are damped for non-zero Lewis number, it is reasonable to 
keep only those modes which are driven directly by the convective amplitude. 

For free-slip-permeable boundary conditions it turns out that only a single mode C is 
excited to lowest order. The form of the evolution equation for A and C can then be derived 
without any detailed calculation and irrespective of the boundary conditions by using the 
translation and rehection symmetries of the system. One obtains to third order 


dtA +{s + S 2 C)d^A = ddlA + {a + fC + + hd^C)A + + ..., 


( 1 ) 


dtC = 5dlC -aC + + (hi + h^C)\A\^ + ihi{A*d^A - Ad^A*) + .... (2) 


For these boundary conditions the coefficients have been calculated explicitly from the 
Navier-Stokes equations [Q. Due to the purely advective nature of the nonlinearity in 
the Navier-Stokes equation, no nonlinear terms of the form G” arise in (|^,|^). Since G is 
even under rehections in x no terms like cIrG"' can arise. Note that in l37|J3^,|40ll the term 
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proportional to /i 4 has been given incorrectly. It describes the wave-number dependence of 
the mixing term hi\A\^. A somewhat generalized derivation of is presented in appendix 

For realistic boundary conditions more modes may arise. Considering a complete set of 
vertical eigenfunctions it is in general to be expected that a whole range of them 

is driven simultaneously. However, their damping due to vertical diffusion (cf. coefficient a 
in (^) grows quadratically with their number n of nodes. In addition, their effect on the 
convective amplitude depends on the strength of the relevant coupling coefficients (e.g. / 
in (|I])). Since the coupling coefficients are proportional to the projection of the respective 
concentration mode on the convective amplitude they decrease rapidly with the number of 
nodes. Thus, high-order modes will have only a weak effect. On that account the most 
relevant modes will be 0^°^ and In addition, the modes excited to leading order will be 
odd in the vertical coordinate, since the convective mode is even and the nonlinear operator 
is odd (cf. ([Al|)). Thus, although the mean concentration mode experiences no damping 
at all (it is actually conserved) it may still be less relevant than since the latter is excited 
at a lower order. 

The above arguments suggest to study extended Ginzburg-Landau equations for a con¬ 
vective amplitude A coupled to a single, weakly damped, odd concentration mode = C. 
The importance of such an antisymmetric concentration mode has also become apparent in 
full numerical simulations with realistic boundary conditions Since C is invariant 

under translations and reflections in the horizontal coordinate it can couple to all terms in 
the equation for A as demonstrated in (P to third order. It is expected, however, that it 
will affect A most strongly through the change of the small linear growth rate, which is one 
of the expansion parameters in the problem, and the linear frequency via the term fCA. A 
striking effect of the concentration held on the local buoyancy has been identihed in p^ , pT| . 
In the following I will therefore consider the somewhat simplihed equations 


dtA + sdxA = ddlA -I- (a -I- fC)A -f c|ApA -|- p|A|^A -|- giAd^lA]"^ -f g2\A\‘^dxA, 


( 3 ) 
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dtC = 5dlC -aC + hi\A\^ + Mxl^r. 


(4) 


This version seems to capture the essential mechanisms of the experimental system. To 
account for the subcritical nature of the bifurcation hfth-order terms as well as cubic gradient 
terms have been included^. Note that (|^,|) have to be considered as reconstituted equations 
0 ; Tey are obtained by combining the solvability conditions arising at various orders. 


Therefore not all terms are of the same order (cf. and h 2 dx\A\^). 

In the analysis I will focus in particular on the effect of the gradient term h 2 dx\A\^ which 
expresses the advection of the concentration held by the traveling wave. As shown by Barten 
et al. the wave generates a concentration current which is antisymmetric with respect to 
the mid-plane of the convection layer [^,^]. If the convective amplitude varies in space 
such a current increases the concentration in the top half of the layer, say, and decreases 
it in the bottom half, thus generating an odd concentration mode C. The other gradient 
term ihi^{A*dxA — AdxA*c) will, however, be neglected; for free-slip-permeable boundary 
conditions the coupling coefficients hi and ^4 are both of the order of the Lewis number and 
therefore small. Since the term involving h/i expresses the dependence of the mixing term 
hi|y4p on the wave number it represents a correction to a term which is already small. 

Note that the diffusive term is not due to molecular diffusion, which contributes 

only at higher order (cf. (|A25|) ); instead, it arises through the effect of large-scale variations 
in the concentration held on the local buoyancy of the huid which generate vorticity which 
in turn advects the basic (linear) concentration prohle. 

It should be emphasized that the nonlinear coefficients of the equation for A are not 
the same as those in the usual Ginzburg-Landau equation since the latter is obtained from 
(0,^ in the limit of large Lewis number (|a| -C a) by eliminating C adiabatically. This 
leads to additional contributions to the coefficients of (|lD. The elimination is, however, not 


^In 1^] no quintic terms had been displayed although the saturating term had been used 

in the numerical simulations. 
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completely trivial. If one were to take only the terms to the order given in (|^) one would not 
recover the usual values of the coefficients (not even their limit of vanishing Lewis number, 
cf. (^)). This is due to the fact that terms of higher order in the Lewis number, which have 
been dropped in (^), become relevant. In appendix ^ a derivation of is given in which 
those higher-order terms are retained. The adiabatic elimination of C from these equations 
leads then to the usual coefficients in the Ginzburg-Landau equation. 

The extended Ginzburg-Landau equations ([5,§) are not only relevant for binary-mixture 
convection. They describe quite generally the interaction of an unstable short-wave mode 
with a stable long-wave mode. Such an interaction arises, for instance, in oscillatory convec¬ 
tion with a free surface or for waves on the interface of two immiscible fluids in a channel. 
There the interface position represents the long-wave mode. The latter case has been inves¬ 
tigated in detail in The interaction of short- and long-wave modes is also important 
in waves traveling on propagating fronts. There the front position represents the long-wave 
mode. In the presence of translation symmetry its position itself is irrelevant for the dy¬ 
namics; only its gradients can enter the equations. In contrast to the case discussed here, 
the relevant long-wave mode V = d^C transforms then like V —>■ —V under reflections in x 
Q. In spirit, the extended Ginzburg-Landau equations (|^,0) are related to those describing 


patterns in the presence of Galilean invariance or coupled to a weakly damped mean 
flow 


III. THE CONCENTRATION MODE GENERATED BY A PULSE 

To lowest order the equation for the concentration mode is a linear diffusion-advection 
equation with damping and the concentration mode is driven by an inhomogeneous term. 
It is therefore given by the integral over a suitable Green’s function. Unfortunately, this 
integral appears to have no simple closed-form solution. In this section two limiting cases 
are discussed. 

In order to focus on solutions with a steady envelope the extended Ginzburg-Landau 


9 







equations are written in a frame moving with the pulse, 


dtA + (s — V)dxA 


dtC - VdxC 


ddlA + (a + fC)A + c\A\^ A + p\A\^ A + giAdx\A\^ + g 2 \A\^dxA (5) 
5dlC - aC + h^\A\^ + h2dx\A\‘^ (6) 


For steady solutions of (^) the Green’s function G{x,x') is given by 

G{x, x') = , ^ - x') + - a:'))) 

V + AaS ^ ' 

with 


(7) 


A^i,2 = ^ (VG2 + 4a6 ± V) (8) 

and 0(x) denoting the Heaviside function. The concentration mode can therefore be written 
as 


G{x) = f G{x,x') (^hi\A{x)\'^ + h 2 dx\A{x')\'^^ dx' 


(9) 


VH2 + a5 


rOQ 

/ (hi + h 2 k 2 )e~^^^ |y4(a; + x)\^ + (hi — h 2 ki)kie~^^^ \ A{x — x'^'^dx . 
Jo 


For the perturbation calculation in sec. 0 this integral will have to be evaluated for |H(a;)| = 
Asech(Aa;) +h.o.t. Since I was not able to obtain a closed form for that integral I discuss two 
limiting cases. 

For large a the decay rates ki become large and it is sufficient to expand \A[x — a;')p 
around x. One then obtains 


1 / 1 1 


72+1 


= -\A\^ + 

a 


+ A] 41.41= 




a 


^ + -(^ + -)) ^ + 0 ( 4 ). 

a a I a \ a I 


( 10 ) 


Thus, the concentration mode changes the coefficients of the nonlinear gradient terms in 
the equation for the convective amplitude A (or introduces such terms if they should have 
been absent). It is noteworthy that the contributions to these terms depend on the velocity 
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of the pulse. In this limit, however, the concentration mode can follow adiabatically the 
convective amplitude and one obtains essentially the usual Ginzburg-Landau equation. 

Another, more interesting limit is that of large velocity (or equivalently small a and 
small 6). In this case the time scale set by the velocity of the pulse is short compared to the 
relaxation time of the concentration mode. The concentration mode can therefore not follow 
adiabatically. As a consequence the spatial decay rate corresponding to the decay ahead of 
the pulse is large whereas that describing the decay of the concentration mode behind the 
pulse is small. For V > 0 the rates are given by 

= y b + ^ + 0(1^-“)^ , + (11) 

The slow decay behind the pulse allows an expansion of the exponential and one obtains 

^dx'Y (12) 

For 1/ < 0 one has 

k,^-y + 0(V-^), = + ( 13 ) 

and 


poo 

(hi + h 2 fc 2 ) / x'^\A{x + x)\ 

n =0 -^0 


{ ^ poo 

{hi - h 2 ki)Y,{-kiT / x'^\A{x - x)\^ dx + 

n =0 -^0 


1 


(hi + h2h)j:jj^lrjA(x)\^\. ( 14 ) 


After inserting the expansions for ki into ([I^) one obtains for the concentration mode 


~ + ^dM{x)\^ + 

hia roo 


\Mx + Sx)\^dx' 


1/2 


poo 

/ x'\A{x + Sx')\‘^ dx'+ 0{—) 
Jo 


(15) 


with S = |1/|/1/. This expansion is not valid uniformly in space since the gradients in C{x) 
become small far behind the soliton without C{x) itself being small. The damping term 
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therefore becomes important and ensnres that C{x) goes to 0 snfficiently far behind the 
pnlse. This can be taken care of by a matching procednre. For the pertnrbation calcnlation 
below this is, however, not necessary since the soliton amplitnde decays exponentially there. 
Since the envelope of the soliton of the nnpertnrbed nonlinear Schrodinger eqnation is even, 
the odd contribntions to C will affect the velocity of the pnlse. In particnlar, for hi = 0 one 
expects a relationship between the gronp velocity s and the pnlse velocity V of the form 

V ^ s + + 0(V-^), (16) 

where the coefficient K depends on the amplitnde of the pnlse. This is conhrmed by the 
pertnrbation analysis below in which K is calculated in detail. Thus, in the presence of the 
concentration mode the pulse velocity is a strongly nonlinear function of the linear group 
velocity s and of the other coefficients of the extended Ginzburg-Landau equation. Eq.fpTl) 
suggests even the possibility of a multi-valued connection between V and s and a hysteretic 
transition between the different solution branches. This rich behavior is to be contrasted 
with that described by the complex Ginzburg-Landau equation alone for which the pulse 
velocity is strictly linear in s even if nonlinear gradient terms are included p^ . 

IV. THE SOLITON EVOLUTION EQUATIONS 
A. Perturbed Nonlinear Schrodinger Eqnation 

In this section the evolution equations for the solitary wave are derived by considering 
it as a perturbed soliton. The nonlinear Schrodinger equation has four continuous symme¬ 
tries: translations in space, phase shifts, scaling of the amplitude and Galileian invariance. 
The soliton solutions therefore form 4-parameter families of solutions characterized by their 
location xq, phase 0o, amplitude A and wave number q. Under small dissipative and dis¬ 
persive perturbations, as they are introduced by the extended Ginzburg-Landau equations 
(^,^, the dynamics is predominantly within one of these families and can be characterized 
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by the slow evolution of xo(t), (/>o(t), X(t) and q{t). The corresponding evolution equa¬ 
tions are obtained from solvability conditions arising in the perturbation expansion. In the 
nonlinear Schrodinger equation these are due not only to proper but also to generalized 
zero-eigenvectors of the linearized operator in question P9|,|^ . 


Thus, I consider the perturbed nonlinear Schrodinger equation 
idtA + i{s - V)d:^A + ]^dlA + \A\‘^A = eP, 


(17) 


where P contains the remaining terms of the extended Ginzburg-Landau equations including 
the contribution from the concentration mode through eq. ( 0 ). Since the explicit expressions 
for the concentration mode were only obtained in a frame in which the soliton is at rest, 
eq. (pT|) has been transformed to a frame moving with velocity V and it is this velocity that 
needs to be determined. The amplitude A is expanded around the single-soliton solution. 


A(x,t) = (/lo(e) + e/li(e,T) + ...)e«l'^>''+‘'», 


(18) 


with 


Ao = A(T)sech(0), 0 = \{T)9 

and 


(19) 


dt» = -v{T), d,» = i, a,<i, = u,{T), a,^ = o. (20) 

The quantities characterizing the soliton within the family are allowed to vary slowly in time 
in order to eliminate the secular terms arising in the perturbation expansion. In addition, 
they - as well as G - are expanded in e. 


A — Ao -l- eAi -|-... 

, q = qo + eqi + 

V = Vo P evi P ..., 

(21) 


U! = UJq P 6UJi 

V = VoPeViP ... 

(22) 


Note that the perturbation P can also change the width of the soliton. Therefore a long¬ 


wave theory, in which 0 is a slowly varying function of x and t is not sufficient [52|. In the 
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present ansatz any time-dependence of A and q leads to diverging time derivatives of the 
phase for large values of 6. This poses, however, no problem here since the (bright) soliton 
goes to 0 exponentially fast there. Inserting the ansatz (p!8D-(p^) into (^) yields at lowest 
order 

fo = s — Vo -|- go, Uq = ^o)- (23) 

Thus, to this order the amplitude and the wave number as well as the velocity and the 
frequency of the pulse are undetermined. At 0(e) one obtains 

idtAi + CAi = idtAi — —AqAi -|- —Aq^IjAi -|- Aq{2Ai + A\) = J, (24) 

where the inhomogeneity I contains in addition to P the contributions from the derivatives 
of the slowly varying quantities. The operator iC is singular. It has two proper zero- 
eigenvectors IAq and ^eAo, which arise from the translation and phase shift symmetry 
of the unperturbed nonlinear Schrodinger equation. In addition, iC has two generalized 
eigenvectors iQA^ and ©c^eAo -|- Aq 

iC {i0Ao} = —AqC?0Ao, iC {©cleAo -I- Aq} = *AqAo, (25) 


which arise from the scale and the Galileian invariance. In general, Ai will contain secular 
terms. Weinstein has shown |^9| that Ai remains bounded for times of 0(l/e) if I is 
orthogonal to the generalized null-space of in which case Ai is also orthogonal to that 
space. To project il onto that null-space an appropriate scalar product is dehned by 


(A, B) = Re A*{e)B{e)de^ . (26) 

It corresponds to the usual scalar product on 7^^ once the complex quantities are separated 
in real and imaginary parts. Using the appropriate left eigenvectors, the four solvability 
conditions are given by 


Re / fA q/ d© = 0, 


(27) 
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/ oo 

OqAqI dQ = 0, 


( 28 ) 


J — OO 


poo 

Re / iOAQl dQ = 0, 


( 29 ) 


poo 

Re {A„ + edeAa)Id0 = O 


( 30 ) 


— OO 


Inserting the expansion ([T8|) one obtains 



0 


( 31 ) 


( 32 ) 


( 33 ) 


e>i = AoAi + (/o(Vi +t'i) 



1 r°° 

sech(0)P^ dQ + -p / © sech(0) tanh(0)Pr dQ 


/ 


( 34 ) 


with P = Pr + iPi- Eqs . (1^ , 13^ yield evolution equations for the lowest-order contributions 
to the amplitude and to the wave number. This determines also the velocity and frequency 


to lowest order (cf. (^)). The quantity of most interest in the present context is the effect of 


the concentration mode on the velocity of the soliton. To lowest order it is determined by the 
group velocity s and the wave number go, ho = -s + go- The relevant question is therefore how 
the perturbations affect the wave number. It turns out that spatial variations in the growth 
rate due to the concentration mode, which were also identihed in the numerical simulations 
of the Navier-Stokes equations, lead to go = 0 independent of their strength. Therefore, no 
change in velocity arises at this order. One then has to determine Vi which in turn depends 
on gi (cf. (^31) ). This contribution to the wave number is, however, not determined until 
O(e^). Spatial variations in the frequency on the other hand affect the wave number already 
at lower order. 

In order to go to O(e^) one needs to determine Ai. This can be done by variation 
of constants using the homogeneous solutions of the real and the imaginary part of 
respectively. 


A^^^ = sech(0) tanh(0). 


(35) 
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= 0 sech(0) tanh(0) + - cosh(0) — sech(0), 


\{h) 
hr 

^50 = sech(0), 

^ 2 ? = 0sech(0) + sinh(0). 

The general solution is then given by 


( 36 ) 

( 37 ) 

( 38 ) 


Ai — (Fir + Glr)Air + (T2r + G2r)^L^ + 
i ({Fli + Gii)A\^ + (F2i + G2i)^L^) 


(39) 


with 


Fir = ^J IrAt^dQ, 

IrAf^ dQ, 


For — 


■'0 

-3 




Fu = 


0 • 
1 

T 2 


liAfde, 


F 2 . = 


( 40 ) 

( 41 ) 

( 42 ) 

( 43 ) 


and I = Ir + ili- The coefficients G 2 allow the elimination of certain exponential divergences 
of y4i for large |0|. The solvability conditions arising from the projection onto the proper 

I eliminate the other exponential divergences. The contributions 


zero-eigenvectors 


from the coefficients Gi lead only to a hxed shift of the position and of the phase of the 
soliton. They can therefore be set to 0. After inserting Ai into ([T7|) the required equations 
for Ai and qi are obtained by applying the solvability conditions at O(e^). 

B. First-Order Evolution Equations: Effect of C on the Frequency 


The extended Ginzburg-Landau equations (0^) lead to the perturbation terms 


P = i(drdlA + arA + clAI^A + p\A\^A + giAd^\A\^ + gMfd^^A + fAG). 


( 44 ) 


Here G is given for a and 6 small as compared to V by 
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c = 


/lo I 4 i9 hoSX ^ . ..9 


(^ + Vx) ir + ^/e““ (® - ®') l•^<®')l'''®' 

with A = |A|/A. In view of other possible applications of this perturbation expansion a 
slightly more general perturbation is used 


P = i{drdlA + ttrA + c|ApA + p|A|"^A + (pi + iai)Adx\A\‘^ + (p2 + i(y2)\A\‘^d^A + 

1 fSAoo 

{po + iao)A— / \A{Q')\'^dQ'+ 

A Je 

1 fSAoo 

(p_i+^a_i)A- / (0-0')|A(0')pd0'. 

A^ Je 


leads to 


Insertion of this perturbation in eqs.( 

d ^ 16 AqP^ /4cr 2 dr 4 ( 72 ^ 0 ' 

dT^° ~ 15 ^ V~3 3 

d 4:dr\lqo 2 A^(To /4(T2 SaA 4 

-go = ^AAoa_,-^-^ + T^j ^0, 


Aq + 2SA XqPo + ^2 — 2 dr — 2 Aq, 


t>i = a - V - + ?|l) Aj - p„ + A 

8A^p, 


P-lT^ 
6 An 


Ui = 


+ (q + p2qo) Aq + (Ai + ASao) Aq + go('^^i + 4h) ~ 


0--1 


(46) 

(47) 

(48) 

(49) 

(50) 


Eqs.(|47| j48|) constitute two coupled nonlinear evolution equations for the amplitude and the 
wave number of the pulse. This could in general lead to interesting dynamics. The quadratic 
term sSAX^pq in (|47D indicates that the subcritical pitchfork bifurcation is perturbed to a 
transcritical bifurcation. 

Here I focus on a discussion of the influence of the concentration mode. On the one hand 
this restricts the choice of the coefficients aj = fiPj, j = —1, ...2. On the other hand, the 
coefficients depend on the velocity V of the pulse (cf. (|T^) which in turn depends on the 
wave number. This introduces additional complexity. Explicitly, one obtains 


Aa = 


ISpr-Ap 

15 


5 ((“t ^ W) ^ 


X3 


+2A^ 



frXl + 2 [ttr- drql - 


frhia 

Vi . 


Aq + 0{Vq ^), 


(51) 
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d 

dT 


qo = - 


4 /2M 

15 I 


fi + 92,1 + 2 Qi^i ) Aq — 


((hi ah2\ \ 2 , >^ofihiaAS 


'0 


(1 ASa7T‘^\ {2 Xld + 3 a) frh 2 


3 Vi 


(2(y'l,r + 92 , r) Aq + 0(Vq ^), 


fXoAS , 2A2^-a^^ J Ag , AqA^ 


fih2 + 


Qoi'Vi + Vi) + 92,rX^qQ + AqAi + 


8 ApPi 

9 


+ o(V)- 


( 52 ) 


(53) 


(54) 


Eqs. (|5l|j5^ ) show that Aq and go and therefore also Vq are 0(l)-qnantities and can vary over 
an 0(l)-range. The smallness of the pertnrbations expresses itself only in the slowness of 
the dynamics. Note, however, that the solvability conditions are only valid in a frame of 
reference in which the pnlse is steady, i.e. n = 0, ^A = 0 and ^g = 0. This limitation 
arises from the expansion for the concentration mode (|T^). 


Since go = Vq — s eq. (|5^ can be viewed as an eqnation for Vq. Solving for the gronp 
velocity s one obtains 


s = l/n 


fjhi /J_ 
2 dr IVq 


3 SAa A 1 

2 Wj ^ 


Vidr 


-a + -ApA ) + 

r-Z\ 


^(2gii + 92i)Xl + OiV^ ^). 


(55) 


Thns, the pertnrbation terms arising from the dynamical eqnation for C lead to a strongly 
nonlinear connection between s and the pnlse velocity Vo (0)- This is to be contrasted with 
the resnlt one would obtain from a complex Ginzburg-Landau equation alone in which the 
coefficients of the nonlinear gradient terms are hxed [Q, 


s = En- 


3 AAct-i 
4 Xodr 


H-^(cT 2 + 2 (Ti). 

2 dr 5d/ ’ 


(56) 


In 


the perturbation terms lead only to a hxed shift in the velocity (as do the terms 


^Note that to leading order in I/Vq the amplitude Aq is independent of s and Vq. 
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involving gj in (^)). Although the shift depends on the amplitude Aq it does not give a 
natural explanation why Vq should be close to 0 over a range of parameters. 


By contrast, both the terms proportional to hi and to /i 2 in (^5]) are suggestive of a range 
of parameters with small Vq. This is indicated in fig.|^ where the contributions to Vo(s) from 
the hi-term (dashed line) and from the h 2 -term (solid line) are sketched separately for 
fihi > 0 and fih 2 > 0 [^. Of course, the singularities at Vb = 0 are unphysical and arise 
from a break-down of the expansion (|T^ for small Vq. Strictly speaking, even the saddle- 
node bifurcations suggested by fig.Q are beyond the present expansion. In that regime the 
diffusion and the damping term in (|^) come into play. It is expected that they will lead to 
a smooth connection of the two branches of Vb(s) and thus to a whole range of parameters 
(with s > 0) in which Vq is small. This is confirmed by the numerical simulations discussed 
in Sec.0 below. An analytical description that removes the singularity would require a better 
approximation of the concentration mode than that given by (^). Obviously, for negative 
values of fihi and fih 2 the concentration mode leads to a speed-up of the pulse (for s > 0). 

For free-slip-permeable boundary conditions it is found that fih 2 > 0 and fihi > 0 
(cf. ( |A22| , |A25D ). Thus, within that approximation the extended Ginzburg-Landau equa¬ 
tions yield a slow-down of the pulse, as is observed in experiments. For realistic boundary 
conditions the coefficients in eqs.(|^J^) are likely to be different]^. However, to leading order 
in the expansion in 1/Vb only the sign of the two products f^hi and /j/i 2 is relevant to deter¬ 
mine whether the pulse is slowed down. The other coefficients determine only the strength 
of the effect and whether the pulse velocity shows hysteresis. Thus, the mechanism is quite 
robust and is likely to be relevant in the experimental system. It arises from the effect of the 
concentration mode on the local frequency of the wave through the term fiCA. The spatial 
variation of the concentration mode implies a spatial variation of the frequency which leads 


^Note that the nonlinear coefficients in (^) differ from those of the complex Ginzburg-Landan 
equation derived in (cf. appendix]^. 
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to a differential phase winding. Since the gradient term h 2 dx\A\^ as well as the advective 
term Vd^C induce an odd contribution to C, the phase winding results in a change of the 
average wave number go and as a consequence in a change in velocity. 

Through the wave number the perturbation induces also a frequency shift ujq = Ao/2 + 
^ 0/2 — go no (as compared to the Hopf frequency) of the pulse in the lab frame. It is worth 
noting that this change in frequency differs from the frequency shift ootw = -^o ~ 9'o/2 ~ 
of an extended wave at the same wave number. Whether this difference is relevant for the 


fact that the pulse frequencies in the experiments and in the numerical simulations 
are larger than those of the extended traveling waves is not clear at this point. 


C. Second-Order Evolution Equations: Effect of C on the Growth Rate 


In earlier work Barten et al. had pointed out that the experimentally observed slowing- 
down of the pulse could be due to the fact that the concentration held leads to a reduction of 
the local growth rate of the convective amplitude ahead of the pulse [pl|-p^. This observation 
motivated the derivation of the extended Ginzburg-Landau equations (|^^) |^. Eq.(p5D 
shows, however, that changes in the local growth rate do not affect the velocity at leading 


order. Variations in the growth rate come in only at 0(e) (cf. (|53|)). At that order any 
change in the wave number gi becomes relevant as well. Since gi is determined only at 
O(e^) this implies that one has to go to next order in the perturbation expansion in order 
to capture this effect. 

To simplify the calculation it is now assumed that the leading-order, dispersive effect 
on the wave number vanishes, i.e. aj = 0, j = —1,..2. In addition, since hi is small for 
free-slip-permeable boundary conditions it is neglected as well (thus, p_i = 0). Once the 
solvability conditions at 0(e) are met, the amplitude Ai is determined using the general 
solution (pQl). Inserting it into the solvability conditions one then obtains at O(e^) 

d , /16 


PrAg + 2{2cr — dr)XQ + dAS'poAo + 2(a,. — drq^^ Ai 
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X 2176,7 
-AdrqoXoqi + \PrPi + 


/ 2 \ 4 

Aq Up2q + 0 + Ci)pr + -{Cr - dr)pi] + -Aq {{Cr - dr){qoP 2 + q)) , (57) 


32 
15 

d 4,2 J 8 ^ , , 32 , 1 ,,g 

— ~-^\qidr — 2«rAlAogo + -^Pr{pl + gP2)Ao 


/4,8 4,1 128 ^ 136 \ ,4 

+ ^ 15^3^^ ~ ~ -^QoPidr + -^PoPrj Aq 

4 2/1 

'^^oPoPl-^S + - f -(4Cr + dr)po + 2 (Oj, — drql) — 4 {drCiqo + P2<^r?o) ) A, 
+2AoAS'pq + 2{ar — drqQ)po. 


(58) 


the equations for V 2 and uj 2 contain A 2 and q 2 which are only determined 


Similar to eqs.(M 
at 0{e^). 

For the equations describing the effect of the concentration mode the coefficients pj have 
to be inserted. Eq.(^) shows that without dispersion qo = 0 and as a consequence Vq = s. 


For s = 0(1) the equation determining vi will be linear in vi (cf. (|49|) ) and will therefore 
not be able to capture the hysteretic behavior found in previous numerical simulations . 


I therefore take s = esi. In order for the concentration held to remain a small perturbation 
one then has to take ^2 = 0(e) as well. In addition, for the expansion (plSf ) to be valid a and 
6 have to be small as compared to 17 = eVi. With these changes the solvability conditions 
at 0(e) yield again (pl| ) -(|5^) with Vi replaced by Vi — si and Vq replaced by Vi. At O(e^) 
one obtains 


d , (. „ , ( Aq A5aAo\ , Idp^Af 


dT'Ai \ Vi 


+ 


-|- 2 (2 C). — d-p) Ag -|- 2 (flj. — d^q^) j Ai 


2176 XlprPi 64 


-4d^gogiAo + + ^{cr - dr)piXl 

945 15 

QAfrPih2Xl AfrV2h2Xl 


45 lA 


+ 


3 I/A 


+ 0(Ef^), 


(59) 


d 4 , ,2 ^ j \ \ 

^91 = “gWrAgl?! — -dpXiXoqo 


128 


qoPidpXQ 


(2 ,^^2 7 2\ 136prAn\ h2/r« 

+ f g (4 Cr + dr) Ag + 2 (Or - dr^g) H-- 


+ ( o (or — dr^o) Ag + 


32prA^ 4 
63 ^ 15 V3 


Cf df I Q 


A h2frd 




+ 0(^4 


-3^ 


( 60 ) 


In the steady case, one obtains from (|^ , |5^) go = 0 and 
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SPrAo * , , „ , 

Ur o I 2, Cr 

15 ^ V ^1 


2/r^2\ ,2 fr^Sah2Xo , 


A^- 




+ 0(n- 


(61) 


This yields for Ai and qi 


Ai = 


16 PjAq (34p^Aq — 21 dr + 21 Cr) 


63 16prAo — 5 dr + 10 Cr 
16 frPih2Xl (4prAo — 105 dr) 1 


Qi = 


63 (16 PrAo — 5 dr + 10 Cr)^ Vi 

frh2 /26 ,o 1 


+ 0(1^1 ), 


(62) 


r2 3^*"^ 2dr)^ o. 

105d 1/2^ "^0 (SPr-Ao + 7(Cr - dr)) 6 + 0{V^ ^). 


drV^ 

frh2 


To determine the velocity Vi the expression for qi is inserted into 
Using (l6l|) to replace Aq one obtains 


si = Ui 


26 Pr Aq 1 a H ri I f U A ^2 2 O, 

IbdT + 3 0 + zj r l35 0^ “ zj ■ 7l 


(63) 


with hi = 0 = Cl- 


h2fr 


J Vi 

+0{V,-^) 


(64) 


Since the amplitnde Aq appears only in the terms of order 0{Vi^) it can be replaced by the 
solntion of (^) in the limit Vi —>■ cx). In that limit (|5TD describes an nnpertnrbed backward 
pitch-fork bifnrcation. Of particular interest is the amplitude at the saddle-node bifurcation 
and at the Hopf bifurcation (or = 0), 

5(dr — 2 Cr) 5/rh2 


A 


0,SN 


+ 


+ 0(1^1 ), A^,^ = 2A, 


2 

0,SN- 


16 Pr SprUi 

At these special points the velocity is related to the group velocity si by 


( 66 ) 


Si — Vi,5Ar -1- 
Sl = Vi^H — ( 


(4 Cr -|- 61 dr) I (14 Cr -|- 53 dr) ^ h2fr , 


^ 336 (—Pr) 

(2 Cr — dr) 


(2 Cr — dr) d -|- 


120 ^ ■''"2 


drUi^ 


+ 0{V,%), (66) 


SN 


28 (-Pr) 


(Cr — 11 dr) 5 -\- 


(2 Cr — 11 dr) a\ h2fi 
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(67) 


To interpret (^,^,^) note that pr < 0 and 2cr — dr > 0 for the pulse to exist stably in the 
absence of the concentration mode, / = 0 (cf. (0)). Therefore, if frh2 > 0 the concentra¬ 
tion mode slows the pulse down when it is first created in the saddle-node bifurcation. This 
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condition, which is satisfied by the coefficients in the free-slip-permeable case [Q, corre¬ 
sponds to the concentration field being advected in such a way that the contribution to the 
buoyancy from the concentration is reduced ahead of the pulse. Eq.(^) supports therefore 
the physical interpretation for the slow-down given previously 0-0.0- 

With increasing the situation can, however, change. Eq.(^) shows that the contribu¬ 
tion to the slow-down from the damping term proportional to a decreases with increasing 
amplitude Aq. The contribution from the diffusive term proportional to 5 can have either 
sign; while the term proportional to always leads to an increase in the velocity with 
increasing a^, the effect of the term proportional to Aq depends on the ratio d^/cr- For 


Cr > lldr 


( 68 ) 


its contribution increases the velocity, whereas in the opposite case the ‘trapping’ is enhanced 
with increasing amplitude. In the immediate vicinity of the saddle-node bifurcation the 
amplitude Aq increases more rapidly than and will dominate the behavior. Thus the 
outcome depends on the relative strength of damping and diffusion; the pulse velocity will 
increase if 


a 15 dr Cr\ 

6 91 {—Pr) V dr) 

and decrease otherwise. Thus, if > lld^ the pulse velocity increases with independent 
of the strength of the damping term and, in fact, for sufficiently large the concentration 
mode even leads to an over-all acceleration of the pulse (cf. (^) at = 0), i.e. the pulse 
is even faster than the linear group velocity. If (|69|) is satisfied but not (|68|) , i.e. Cr < lldr, 
then the velocity still increases with but remains below the group velocity as long as 

ttr < 0 . 



These results show that the response of the pulse to the concentration mode is more 
complex than suggested by the original picture in which the concentration mode only led 
to a ‘barrier’ of reduced buoyancy which can ‘trap’ the pulse Since the velocity 

of the solitary pulse is closely connected to its wave number q (cf. any effect of 
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the concentration mode on the local wave number becomes just as important. In fact, 
the dependence of the local frequency on the concentration mode, which changes the local 
wave number, selects a velocity of 0(1) (cf. (^)). But even when this dispersive effect 
is absent (/* = 0) the variation of the concentration mode can select a non-zero wave 
number [^| somewhat similar to the selection of the wave number by ramps in the control 


parameter M,|^. Eq.(|^ shows that in the absence of any change in wave number the 
barrier generated by the advected concentration mode, which is represented by the term 
proportional to ^ 2 , would always decrease the pulse velocity. 


V. NUMERICAL SIMULATIONS 

The perturbation analysis in sec. [IV B| and sec. |IV Cj revealed two distinct mechanisms 
that can lead to a slow-down of the pulse due to the coupling to the concentration mode. A 
dispersive mechanism which affects the velocity through the wave number and a mechanism 
that is based on a suppression of convection ahead of the pulse. The analysis suggests that 
there is a range of parameters over which the pulse is slow. However, due to the expansion 
employed for the calculation of the concentration mode the analysis is strictly valid only for 
large velocities and can only be suggestive with regard to the behavior at small velocities. 
To go beyond this limitation the extended Ginzburg-Landau equations ii) are solved 
numerically. 

To study the dispersive ‘trapping’, eqs.(^,§) are solved with / purely imaginary. The 
values of the coefficients are chosen as dr = 0.15 + i, a = —0.24, / = 0.5q c = 2.4 -|- 2i, 
p = —1.65 + 2i, 6 = 0.1, a = 0.08 and the boundary conditions are periodic. Fig.^ shows the 
velocity of the pulse for hi = 0, h 2 = 0.3 and for hi = 0.5, /i 2 = 0. As expected from (^) the 
pulse is slowed down by the concentration mode (the dotted line gives the velocity without 
concentration mode, U = s). The simulation bears out even the hysteretic transition to a 
‘trapped’ pulse which is suggested by (|55[ ) but outside its range of validity. This transition 
can be understood intuitively by noting that the ‘amount of G’ generated by the pulse at 
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any given location depends on its velocity. If the pulse is slow C{x) becomes large and slows 
it down even more. If the pulse is fast, however, C{x) is small and the effect on the velocity 

is small. 

The terms hijAp and h2dx\A\^ differ in their symmetry and their physical interpretation. 
Whereas the h 2 -term corresponds to the advection of the concentration mode by the wave, 
the hi-term represents the generation of C. For the concentration mode to affect the (over¬ 
all) wave number it must have an asymmetric contribution. For n = 0 the hi-term leads to 
a symmetric C-prohle. Consequently the wave number is increased on one side of the pulse, 
but decreased by the same amount on the other side. The over-all wave number and the 
pulse velocity remain unchanged. For non-vanishing velocity, however, the concentration 
held lags behind and changes the wave number as illustrated in hgj^. Due to its symmetry 
this term can induce ‘trapping’ of the pulse independent of its direction of propagation; for 
/i 2 = 0 (^) is invariant under s —> — s and Vb —>■ — Vq. 

By contrast, the h 2 -term breaks the rehection symmetry at any velocity. For fih 2 > 0, the 
concentration held pushes the pulse to the left independent of its direction of propagation. 
Thus, for positive group velocity s the pulses are slowed down, whereas for negative s they 
are accelerated. 

If the dispersive ehects are small the pulse velocity is ahected at next order through the 
change in the local growth rate given by frCA. This is the mechanism which was pointed 
out by Barten et al. [pT| -p^ . Some numerical simulations of the extended Ginzburg-Landau 
equations (|^,^) for this case have been presented previously |l37|j5^ . As suggested by (^) the 
qualitative dependence of the pulse velocity on the group velocity is similar to that due to 
dispersion. The results of simulations for a = —0.24, d = 0.15-1-i, c = 2.4-|-2i, p = —1.65-1-2i, 
6 = 0.1, a = 0.02, / = 0.5, /i 2 = 0.3 are shown in hg.|. Again, for sufficiently large group 
velocity a jump-transition from fast to slow pulses is observed. In the vicinity of the jump 
the trapped pulse becomes unstable to oscillations. It becomes more extended - in particular 
at its trailing end, where the concentration mode leads to an increase in the local growth 
rate - and phase slips arise. This is shown in hg.|. Here a pulse is depicted at the time 
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of the phase-slip process which is marked by the localized amplitude depression and the 
diverging local wave number. The phase slips lead to an oscillatory behavior of the velocity, 
the extrema of which are indicated by the dashed lines in hg.^. When the group velocity is 
increased further the trailing end develops larger peaks which then grow by themselves and 
eventually convection spreads throughout the system. 

Experimentally, the only parameter that can be varied easily is the Rayleigh number, 
which enters the growth rate ar- An interesting question is therefore whether the ex¬ 
tended Ginzburg-Landau equations allow a hysteretic transition to fast pulses by changing 
ttr- Eq.(|64D suggests this possibility, but is not valid for small velocities. This question is 
investigated for three cases. 

The growth rate enters (Q) in three places. The explicit dependence leads to an increase 
in velocity for all frh 2 > 0. So does the contribution (via the amplitude Aq) from the damping 
term proportional to a. The diffusive term proportional to 6 can, however, lead to a further 
decrease in the velocity with increasing amplitude. If inequalities (|68| , |69D are satished, the 
pulse velocity is predicted to increase in the immediate vicinity of the saddle-node bifurcation 
and reach a value larger than s for = 0. The result of a numerical simulation of this case 
is presented in hg.^ for s = 0.3, d = 0.05 -|- 0.5q / = 0.25, c=2.4 + i, p = —1.65, a = 0.03, 
S = 0.1, h 2 = 0.3. It shows indeed a hysteretic transition from a trapped pulse to an 
accelerated pulse. Interestingly, the trapped pulse becomes unstable to oscillations in the 
vicinity of the transition. Here it is the variation in the peak height of the concentration 
mode which leads to the oscillations in the velocity. Their minimal and maximal values are 
denoted by dotted lines. 

Fig.^ shows the result of simulations for a case in which Cr < lldr but inequality (|6^) is 
still satished (solid squares, dr = 0.3, other parameters as in hg.|^). As expected from the 
perturbation result the pulse velocity increases in the vicinity of the saddle-node bifurcation 
(or ~ —0.5), but it remains below the group velocity up to = —0.15. For larger values 
of ttr the pulse becomes unstable. Finally, if both inequalities (|5^ , |6D|) are violated the 
perturbation analysis suggests that the pulse is slowed down even further with increasing a^. 
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Two such cases are shown in hg.|^ {dr = 0.6). For fr = 0.35 (open circles) and for fr = 0.25 
(open diamonds) the pulse is ‘trapped’ and its velocity decreases up to ar = —0.43. For 
larger values of the pulse starts to grow and eventually hlls the whole cell. For weaker 
coupling to the concentration mode {fr = 0.2, open triangles) the pulse velocity decreases 
with increasing in the immediate vicinity of the saddle-node bifurcation, but for larger 
Qr the velocity starts to increase. 


VI. CONCLUSION 


In this paper I have investigated the influence of a long-wavelength mode on the dynamics 
of a traveling-wave pulse in the perturbed nonlinear Schrodinger equation. The study was 
motivated by experimental and numerical results on convection in binary liquids. It showed 
analytically that such a mode strongly affects the propagation velocity of such a pulse. In 
particular, two distinct mechanisms were identihed which can lead to a ‘trapping’ of the 


pulse. As expected from numerical simulations of the full Navier-Stokes equations this 
can be due to a reduction of the local growth rate of the traveling-wave mode. This effect 
arises, however, only at second order. Already at hrst order the frequency of the waves is 
affected by the long-wave mode which induces a change in the wave number and connected 
with it in the velocity of the pulse. This effect was neither anticipated by the numerical 
simulations nor by the experiments. For simplihed boundary conditions both effects are 
found to reduce the velocity unless dispersion is very large. Within the framework presented 
here hysteretic transitions can occur between ‘trapped’ pulses and fast pulses when the local 
growth rate (i.e. the Rayleigh number) is increased. So far no such transitions have been 
found in the experiments or the numerical simulations of the Navier-Stokes equations. 

An increase in the Rayleigh number in the experiments often leads to long pulses which 
cannot be regarded as perturbed solitons anymore. Instead they can be considered as a 
bound state of two fronts. The dynamics of such fronts is also strongly affected by the long¬ 
wave mode. In particular, within the extended Ginzburg-Landau equations discussed here 
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it has been found that their interaction can change from attractive to repulsive when the 
direction of propagation of the pulse changes sign SB- Thus, backward traveling pulses 
are stable and pulses of different sizes can coexist in agreement with experimental results 
35|. Further studies show that the coexistence of stable pulses found in the numerical 


simulations [^| can also be understood within this framework. 

No attempt has been made so far to investigate the effect of increasing the Lewis number 
(i.e. the mass diffusion) using full numerical simulations of the Navier-Stokes equations. 
According to the present analysis the pulse velocity should increase drastically when the 
Lewis number is increased. Such an investigation would provide a stringent test for the 
relevance of the extended Ginzburg-Landau equations to binary-mixture convection. 

There have been attempts to find stable localized pulses in two-dimensional convection. 


To date these pulses have not proved to be stable but decay after long transients . The 
origin of this instability is not understood. In addition to the long-wave concentration mode 
a long-wave vorticity mode (mean flow) is expected to be relevant. A related interesting 
problem is the dynamics of two-dimensional propagating fronts. Since the concentration 
mode slows down certain types of fronts one may expect that it could lead to a transverse 
instability of such fronts similar to the Mullins-Sekerka instability in solidification . 

The relevance of the equations discussed in this paper is not conhned to convection in 
binary liquids. Equations of the same form have also been derived for waves on the interface 


between two immiscible liquids flowing in a channel |^. It is conceivable that nonlinear 
gradient terms arising in the description of other wave systems are also due to the elimination 
of a second mode (e.g. [^,^]). In regimes in which such a mode becomes slow as compared 
to the wave it would have to be described by a separate dynamical equation and could have 
similarly strong effects on the waves as found in the present paper. 
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APPENDIX A: DERIVATION OF THE EXTENDED GINZBURG-LANDAU 


EQUATIONS 


In this appendix I rederive the extended Ginzbnrg-Landau equations (ECGL) under 
slightly more general conditions than previously |^|. In the previous derivation the Lewis 
number L and the growth rate a of the convective mode were taken to be of the same 
order. In the limit of very small growth rate, i.e. for |a| -C L, the concentration mode 
C can be eliminated adiabatically from the resulting equations to yield the usual complex 
Ginzburg-Landa equation (GGL). However, with the coefficients calculated in [^] one does 
not recover the correct values of the coefficients in the GGL, not even their values in the 
limit of small Lewis number. The generalization presented here overcomes this problem and 
agrees for |a| -C L quantitatively with the GGL. At the same time it shows that the terms 
|Ap and i {A*dxA — Ad^A*), which in the previous derivation were found to have vanishing 
coefficients at quadratic and third order, respectively, become relevant at higher order. 

A simple analysis of (|l|J|) shows that the coefficients derived in do not lead to the 
correct coefficients of the single GGL if the concentration mode is eliminated adiabatically 
for |a| -C L oc a. This is due to the fact that any L-dependence of the nonlinear coefficients 
is shifted to higher-order terms due to the initial assumption L = O(e^) = 0(a). Since 
in the adiabatic elimination the contributions from the nonlinear terms are divided by L, 
higher-order terms of the form L|Ap become relevant. The goal is to obtain certain of these 
contributions without performing the full expansion to the corresponding high order. 

As pointed out in sec.|T|, (|l|J^) are reconstituted equations [@; i.e. the solvability condi¬ 
tions arising at 0(e), O(e^) and O(e^) are not satished separately but are combined into a 
single solvability condition. One could imagine such a reconstitution involving all orders in e. 
The higher-order terms would then have the form L'^A^C^ and suitable spatial derivatives 
thereof. They could be summed up into a term A^C^ (and the respective spatial derivatives) 
with an L-dependent coefficient. Alternatively, one could obtain this L-dependence directly 
by avoiding the expansion in L, thus providing a summation of certain contributions in L 
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to all orders. This is the approach taken here. A more detailed analysis is then required to 
identify for which of the coefficients the procedure restores the complete L-dependence. 


Within the Boussinesq-approximation the basic equations for two-dimensional convection 
in a binary mixture at inhnite Prandtl number read p^|33 


-A2 

-{l+‘ip)d^ 

-09a; 

-Rd^ 

dt-A 

0 

-Rd^ 

dt 

9*-LA 




e 

W 


0 


\ 


(9209a; - d^(j)dz)e 

(9209a; - d^(j)d^){r] + e) J 


(Al) 


where 0 denotes the stream function, 6 the deviation of the temperature from the conductive 
state and r] = {c/'ip — 9) with c being the deviation of the concentration held from the 
conductive state. The separation ratio is denoted by 0 and IZ and L are the Rayleigh 
number and the Lewis number, respectively. For free-slip-permeable boundary conditions 


the helds are expanded as 


0 = e sin( TT ^ ) A + e2{sin( 2ti z) + ^ 

sin(7rz) (^^F^ + iGi + ) e'"™*) + sin( 3 vrz ) 

-bc.c.+ 0(6^), (A2) 

9 = 6 sin( TTz ) AO + e2{sin( 2 vrz ) (^^De + Eg ) 

+ sin( TT ^ ) Qpe + CiiGi - G 2 + Gs) ) + sin( 3TTz)He 

-|-c.c. -|- O(e^), (A3) 

rj = e (e'^"-'‘"'‘*sin(7rz)A0+ Gsin(2 7rz)) +e2{sin(2 7rz) + 

sin( TT z + C 2 iGi + Gs - Gs) + sin( 3ttz)H^ ) 

-fc.c.-I-O(e^). (A4) 


The Rayleigh number is expanded as 

1 _2A i' 1 r 'i 

n = Ro{q) + eR,+e^R 2 with Ro = ^ A (A5) 

At this point the expansion wave number q is performed is left undetermined. This allows a 
consistency check of the resulting coefficients since certain higher-order coefficients have to 
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correspond to derivatives of lower-order coefficients with respect to the wave number q (see 
appendix ^). Eventually q will be chosen to be the critical wave number q^ = 7i/y/2. The 
separation ratio -0 and the Hopf frequency are expressed as 


'Ip = 


ujh = uj{q^ + TT^) with uj = + LVt'^. 


f]2 + l’ 

The linear eigenvector (1, (^ 2 )* of the convective mode is given by 

z(fi2 + l)(l + L)(g2 + 7r2)2 ^ ^ 1 


(A6) 


Ci = 


C2 — Cr 


(A7) 


(1 — iijj)q — iuj 

As mentioned above, the Lewis number L is not expanded. Thus, already at 0(e) a damp¬ 
ing term —dvr^LO arises formally in the solvability condition for C. At this order no time- 
derivative arises which would balance this term. Since L is eventually taken to be small 
(0(e)) it is reasonable to defer the balancing until all solvability conditions for O are recon¬ 
stituted into a single equation. 

At O(e^) one obtains 


Fe 

De 

G 2 


0 ) Frf, — 0 , — 0 , Eg — 0 , — 0 , Ejj — 0 , 

1 (g^ -F TT^)^ (ff^ -F 1) |Ap ^ _ 1 <^) 

2 {L-n‘^-l)Tr 
1 q^ {iL — uj) A R 


16 (112 _|_ ’ 

1 iVF Ti C A 


(AS) 

(A9) 


0,9 — 


jL _ ij \1,±J — UJj _^ li q ^ 

4 (1 -I- L ) o; (g^ -I- 7r2)^ (122 _j_ 4 (-^2 _j_ (-^2 _|_ ^2^3 

1 L + 2 L q^ - i q^ Oj) 

2 a;g (g2 + 7r2) 

1 iq^ ARi 1 7rg^(a;-|-lL)OA 
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1 (—vr^ — n‘^L + 2q^ + 2Lq^ — ?)iq^oj + iuj vr^) A^ 


(AlO) 


H^ = 4: 
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2 a; g (g2 + 7r2) 

i(19i-|-3a;)12^AO 
(122 _|_ ^3 ’ 

i^/2AC( 27 L - 6832 + 1083 ia;) 


Ho = -27 
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(A12) 

(A13) 
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with 
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fci = 513(L^ - - Lif) +3249aJ^ +20496iaj(L+ 1) - 129808 L. 


(A15) 


Note that for simplicity the expressions Hq and have been given only for q = q^. 
In the expressions for G 2 and G 3 , however, the ^-dependence has been retained for nse in 
appendix The amplitndes Gi and Dn can be set to 0 since they only renormalize the 
amplitndes A and C, respectively. Combining the solvability conditions for A and G arising 
at 0(e), O(e^) and O(e^) one obtains the ECGL 


diA -|- (s -|- S 2 G)dxA — dd'^A -|- (a -I- fG -h / 20 ^ - 1 - f^dxG)A -T (A16) 

dtG = 5dlG -aG + hi\A\^ + + h^G\A\^ + 

ihi{A*dxA- AdxA*). (A17) 


The coefficients are given by 
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1 e {00 + iL — i R 2 
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-e TT 


574389 + 28147427+ 51487378 + 25526060 

112 (-^2 + 1)2 ^^3249122 + 182329) 


L + 0{L^), 


h4=-rr 


3 e (112 + 1) tt^ L v/2 


(A27) 

(A28) 


2 112 (L - 112 - 1) 

In these expressions q has been set equal to Qc- Note that hi and are then proportional 
to the Lewis number and are therefore of higher order in agreement with the previous 
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Without a further discussion it is not clear to which order in L the coefficients derived 
above are actually correct; in general the reconstitution does not yield all contributions. 
Consider the expansion with L taken explicitly to be 0(e), i.e. L = eLi. The procedure in¬ 
volves essentially three steps. First the equations are expanded in a straightforward way. For 
the relevant modes A and C this yields evolution equations containing also the other, stable 
Fourier modes like G* etc. In a second step the stable modes are eliminated adiabatically 
using equations like (|A8D-(|A14|)^. This leads to equations of the form 


dtA = A.F(|A|2, G, d,, dt, Li), dtC = ^(|A|2, G, 4, dt, h), 


(A29) 


where the form of the equation is dictated by symmetry. The functions T and Q stand for 
invariant polynomials in their arguments. The notation indicates that A and A* appear only 
with equal powers as required by translation symmetry. The dependence of T and Q on Ri 
is not indicated. Many of the high-order terms differ from the lower-order terms only in their 
powers in Li. By not expanding in L all these terms are summed up into an L-dependent 
coefficient of the lowest-order term with the same structure regarding the amplitudes A and 
G and their derivatives. 

In a third step the time-derivatives in the nonlinear terms are replaced recursively using 


( |A29|) . If only amplitudes and derivatives are counted, this introduces terms of lower order 
since the evolution equations have the formQ 


A ignore in this discussion the possible effect of other critical concentration modes. 

®Of course, the coefficients a and a oa L are small; therefore the true order of the respective terms 
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dtA = aA+ dtC = —aC + ... 


(A30) 


As a consequence, the replacement of dtA contributes to the i?i-dependence (a oc i?i) and 
the replacement of dtC to the L-dependence of the corresponding lower-order term. Even if 
no expansion in L is performed, the latter contributions are not obtained unless the initial 
expansion is taken all the way to the order of the term involving the time-derivative. Thus, 
due to the third step the above procedure will not sum up the L-dependence to all orders 
unless the coefficient in question has no contribution from higher-order terms containing 
time-derivatives dtC. 

Since the elimination of dtC in a given term introduces a factor C, which cannot be 
eliminated by further substitutions, low-order coefficients of terms containing only factors 
of A and its derivatives will not be affected by the replacement of dtC in the higher-order 
terms and are expected to be correct to all orders in L. Coefficients of terms containing a 
factor C, however, may be affected and may therefore only be correct to the order of the 
initial expansion. The terms omitted in (|A3CI|) do not affect this argument since they do 
not decrease the order in terms of A, C and their derivatives. For instance, they lead to 
replacements like dtA —> —sd^A or dtC —>■ hi\A\^. 

The above arguments suggest that the coefficients s, c, hi, ^ 2 , and as given in (|A18|) - 
( |A28| ) are correct to all orders in L. It is therefore expected that the adiabatic elimination 
of C for |a| <C L will yield the same value for the cubic coefficient of A as found in the direct 
derivation of the CGL. 


APPENDIX B: ELIMINATION OF C AND WARD IDENTITIES 

In this appendix it is demonstrated that the coefficients as calculated in appendix ^ are 
consistent with those of the usual CGL in the appropriate limit and satisfy the relevant 
Ward identities H. 


remains the same. 
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To show the consistency with the usnal CGL, C is adiabatically eliminated for e -C T. In 


the previons expansion this did not lead to the correct cnbic coefficient in the CGL since 
for e L higher-order terms of the form LlAp become relevant as well, which are of higher 
order in the previons expansion. They are kept in the expansion presented in appendix 
To wit, hi is of 0{L) and vanished therefore in the previous calculation. In the present 
approach it gives an 0(1) contribution to O, 

9 |A|V (G2 + 1) 


C„, = ^ 


oidMr). 


(Bl) 


8 G2 (L _ _ 1) 

Inserting Cad into (|A16|) gives an additional cubic term in A. The two cubic terms combined 
yield the purely imaginary cubic coefficient of the usual Ginzburg-Landau equation 


Cad 


I TT 

’8^7’ 


(B2) 


in agreement with previous work [61 


An additional test of the coefficients calculated in appendix ^ is possible using the Ward 
identities which relate coefficients of gradient terms with g-derivatives of other, lower-order 
coefficients 10 . These identities arise due to the fact that the Ginzburg-Landau equations 
can be derived by an expansion around any wave number q on the neutral curve. Of course, 
in an expansion around a wave number different than that corresponding to the minimum 
of the neutral curve the group velocity parameter s is complex and the expansion of the 
Rayleigh number contains also the linear term Ri (cf. (^1)). Spatially periodic solutions 
with wave numbers qi = q + Aq differing slightly from q can then be determined in two 
ways. Either q is changed to qi or the difference Aq in wave number is absorbed in the 
slow wave number Q. In the former case the coefficients in the Ginzburg-Landau equation 
change. In the latter case the gradient terms come into play. Any physical quantity P like 
the stream function 0 depends on the expansion wave number q as well as the slow wave 
number Q. Since both approaches have to lead to the same hnal result one obtains the 
invariance condition 


±P(,-eQ,Q)=^. 


(B3) 
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For definiteness, consider a specific Fourier mode of the stream function 0 (cf. (|A2|) ), 

0^^^^ = ^ + Ga + Gg) + 0(e=^)) . (B4) 


From (|A10|,|AlID one finds 


Gj = with gi = const. 


(B5) 


Since Gi can be set to zero one obtains 




(B6) 


with 


P = 


(tt^ - 3g2) 


(B7) 


(g2 ^ 7r2)g' 

The concentration mode G is independent of Q since it is not the amplitude of a periodic 
function. 

Conditions for the coefficients in the extended Ginzburg-Landau equations ( |A16|JAl7[ ) 
arise from taking the Q-derivative of ( |A16|jAl7D and inserting the invariance condition ( |B 6 D . 
It is important to note that the derivative has to be taken at hxed physical parameters like 
the Rayleigh number. Thus, i?i has to be replaced by (Jl — Ro{q))/e. This implies that the 
linear coefficient a{q), which is proportional to Ri, satishes 


- ‘Q) = = C(i)> 


(B8) 


whereas the Q-derivative of all other coefficients is 0(e). Applied to the equation for the 
convective amplitude, (|A16|) , the invariance condition leads then to the expected relations 


. d . d 

s = ie—a,, S2 = t^f. 
dq dq 


(B9) 


The Q-derivative of (|A17|) yields 


d 


2/14 — — (p + P*)hi + 0(e). 

dq 


(BIO) 
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At first sight it might be surprising that is not just proportional to the derivative of hi. 
However, since the eigenvector (1, Ci, which connects the amplitude A with the physical 
quantities, depends explicitly on the wave number q additional terms like (p + p*)hi have 
to be expected in general. 

To check the Ward identities (P^ , |B1U| ) the g-dependence of a, / and hi has to be kept. 
The calculation presented in appendix ^ yields 


a(g) 


/(?) 


hi{q) 


1 i [L — — 1) q‘^ e 1 q^ {uj + i L — iVR) 

8 a;3 (f]2 + 1)2 (^2 + ^2uj{n^ + l) (g2 + ^ 

(Bll) 


1 (27r2L3 + g2^2^2_^4^2^_^2^4^_^2^2^2_^2^4_^2^4) ^^2 

2 0;2 (ff2 + 1) (g2 + ^2)3 + 

. f 1 Ri eVR (fR + 1 — L) L TT q'^ 1 {2 + VR' + q^ VR') q^\ ("gi o') 


1 e i?i g2 (—TT^ L 122 _)_ /^ g2 q 2 _ 2 7^2 _ ^2 ^2 _|_ q,2 ^2^ 

~2 122 ( 1 + L ) (L - 122 - 1) (g2 + 7^2) TT 

TT (122 + 1 ) (g 2 ^ ^ 2^2 ^ 

122 (L - 122 - 1 ) ’ 


(B13) 


where Ri and R 2 have to be replaced by IZ and Rni^q) via IZ = Roi^q) + eRi + e^R 2 . The 
values for s, S 2 and h^ obtained by inserting (pil|,pi^ , |B13|) into (p^ ,|B10|) agree with those 
obtained in (|A19|) and (|A28|) up to terms of order 0(^2) and 0(e) for S 2 and h^, respectively. 


These differences are expected based on the discussion of the correctness of the L-dependence 
of S 2 in appendix ^ and in view of ([B10|) . 
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FIGURES 
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Group Velocity s 

FIG. 1. Sketch of the contributions to the pulse velocity Vq from the hi- and the /i 2 -term (cf. 


(HD). 
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FIG. 2. Numerically determined pulse velocity U as a function of the group velocity s for 
dr = 0.15 + i, a = —0.24, / = 0.5z, c = 2.4 + 2i, p = —1.65 + 2i, 6 = 0.1, a = 0.08 and the indicated 
values of hi and /i 2 . 
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FIG. 3. Typical pulse solution for /i / 0 and hi / 0. The dependence of the frequency on the 


spatially varying asymmetric concentration leads to a shift in the wave number (cf. (52)) which 
changes the velocity of the pulse (cf. (^)). 



Group Velocity s 

FIG. 4. Numerically determined pulse velocity G as a function of the group velocity s for 
a = -0.24, d = 0.15 + i, c = 2.4 + 2i, p = -1.65 + 2i, 5 = 0.1, a = 0.02, / = 0.5, = 0.3. 

The dashed lines indicate the minimal and maximal values of V in the oscillatory regime (solid 
squares). 
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Spatial Coordinate x 

FIG. 5. Pulse solution for s = 0.35 (other parameters as in fig.§). The depression in |^| and 
the diverging wave number indicate that a phase slip has just occurred. These phase slips lead to 
an oscillatory behavior of the velocity. 



FIG. 6. Numerically determined pulse velocity P as a function of the growth rate for s = 0.3, 
d = 0.05 + 0.5i, / = 0.25, c = 2.4 + i, p = —1.65, a = 0.03, 6 = 0.1, /12 = 0.3. The dashed lines 
indicate the minimal and maximal values of V in the oscillatory regime (solid squares). 
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Growth Rate a 

r 

FIG. 7. Numerically determined pulse velocity G as a function of the growth rate for various 
values of dr and fr'- dr = 0.3, fr = 0.25 (solid squares), dr = 0.6, fr = 0.2 (open triangles), dr = 0.6, 
fr = 0.25 (open diamonds), dr = 0.6, fr = 0.35 (open circles). Other parameters as in figj^. Inset 
shows initial decrease in velocity for dr = 0.6, fr = 0.2 as expected from (|^. 


45 














